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ABSTRACT 


The  effect  of  tody  geometry  on  the  equilibrium  turbu¬ 
lent  hypersonic  vake  is  examined  by  computing  the  wakes  behind 
several  bodies  with  different  cone  and  sphere-cone  nose  shapes. 
Typical  re-entry  trajectories  at  several  different  altitudes  are 
used  for  the  free  stream  conditions.  Also,  theoretical  predictions 
for  the  vake  growth  behind  cones  Eire  compared  with  available  bal¬ 
listic  range  measurements. 

In  order  to  determine  the  initial  conditions  for  the 
subsequent  wake  behavior  the  bow  shock  shape,  the  viscous  base 
and  neck  regions,  and  the  inviscid  wake  shock  are  examined .  -  The 
results  for  the  initial  wake  properties  Indicate  that  the  inviscid 
flow  is  made  uoM  by  making  the  body  sufficiently  slender,  even 
thdugh  the  wake  shock  now  gives  a  relatively  significant  contri¬ 
bution  to  the  enthalpy.  Also  it  is ’found  that  in  determining 
the  initial  viscous  core  properties  the  Reynolds  number  is  much 
more  Important  thsin  the  shape  of  the  body. 

From  the  downstream  wake  histories  the  results  indicate 
that  initially  the  slender  body  wake  grows  slower  and  cools 
much  faster  them  the  blimt  body  wake.  Eventually,  sufficiently 
far  downstream  the  wakes  of  all  bodies  should  be  e  unique  function 
of  CjjA.  However,  until  the  asymptotic  downstream  conditions  occur, 
it  is  found  that  even  two  bodies  with  the  same  drag  coefficient 
but  different  shapes  can  have  completely  different  wake  growth  and 
enthalpy  histories.  Rirther,  the  degree  of  this  difference  depends 
considerably  on  the  altitude  (Reynolds  number'.  Since  equilibrium 
electron  histories  follow  the  enthalpy,  a  slender  cone  can  have  a 
very  short  radar  wake  compared  to  a  blunt  body.  On  the  other  hand, 
because  of  the  initially  slower  vake  growth  for  the  cor i  an  even 
longer  time  than  for  the  blunt  body  can  be  required  for  the  radial 
diffusion  of  any  non-reconiblning  foreign  substeince  (e.g.,  ablation 
material)  deposited  in  the  Veike. 
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1.  THTRODUCTIOH 


The  turbulent  vake  behind  a  body  moving  at  hy^rsonie  speeds  can  be 
considered  to  be  composed  of  two  portions  (Fig.  1);  an  outer  inviscid  region 
and  an  inner  turbulent  viscous  core.  For  a  bljint  body  the  large  nearly-nonaal 
portion  of  the  bow  shock  generates  very  hi^  flow  tenq)eratures,  with  the  result 
that  under  equilibrium  conditions  the  outer  inviscid  region  remains  hot  even 
after  isentropic  expansion  to  aaibient  pressure.  Furthermore  the  turbulent  core 
itself}  which  has  undergone  viscous  losses  around  the  body  and  in  the  free 
shear  layer,  can  also  be  hot.  The  subsequent  cooling  process  by  core  growth 
and  turbulent  mixing,  although  much  faster  than  laminar  diffusion,  has  been 
found  (1)^  to  require  distances  of  hundreds  of  body  diameters  downstream, 
largely  because  the  hot  core  is  surrounded  and  fed  by  additional  hot  gas  from 
the  outer  inviscid  region.  The  purpose  of  this  paper  is  to  examine  the  effect 
of  body  geometry  on  the  turbulent  wake.  This  will  be  accomplished  by  computing 
the  wakes  for  bodies  of  vrarious  bluntness  ratios  (nose/base  diameter)  and  com¬ 
paring  with  the  previously  computed  sphere  (l). 

It  is  necessary  to  examine  the  changes  produced  by  body  shape  in 
both  the  outer  inviscid  flow  and  the  initial  turbulent  core  or  neck  properties. 
Since  the  bow  shock  strength  is  nearly  everywhere  quite  weak,  the  inviscid  flow 
changes  considerably  and,  in  fact,  should  be  made  "cold"  by  making  the  body 
slender.  The  contribution  of  the  wake  shock  for  a  blunt  body  was  found  to  be 
negligible  (l)  in  comparison  with  the  bow  shock,  and  mw  must  be  re-examined 
for  a  slender  body. 

Regarding  the  viscous  core  properties  it  is  not  possible  to  state 
a  priori  what  is  the  effect  of  changing  body  shape .  Conventional  methods  of 
calculating  the  boundary  and  free  shear  layers  were  used;  however,  clearly 
more  detailed  analysis  of  the  latter  and,  above  all,  experimental  data  at  the 
higher  velocities  are  needed.  For  example,  the  basic  structure  of  the  wake 
neck  geometry  at  high  velocities  (H  ^  20)  would  be  very  useful. 

Two  Important  physical  processes  determining  the  properties  of  the 
tiu-bulent  wake  are  non-equilibrium  chemlwil  processes  and  body  ablation 
characteristics.  Discussion  and  preliminary  results  regarding  the  former  as 
to  hypersonic  wal»  flews  are  given  by  Lees  (2)  and  Bloom  (3)«  For  the  present. 


Denotes  references  listed  at  end  of  paper. 


we  acBunc  here  chemical  aquilibrluii  of  y\src  air,  except  that  the  specieJL  case  of 
■'pur**''  liffusica  cf  a.  foreign  species  with  zero  rsconbiration  is  caasidered. 

T-'o  iani’.ien  of  todies  art:  ccnsiaeied  here.  These  are  shown  in  Fig.  2 
with  t:ie  pijrtlnc-nfc  fJow  propsrtie-j  given  ta  Tables  I  and  II.  Cse  of  these  families 
consists  of  sphere-ccnes  with  a  constant  cone  angle.  The  second  faaily  are  pure 
cones  of  '."i-ious  angle*.  The  woiwa  are  calculated  for  typical  re-entry  conditions 
of  M  ^  =  20-32  (Table  l)  at  vrorioas  »:titu)k3  (jSeynolde  number  variatlone).  Also, 
thecretici'-l  predi  ctiors  for  the  wake  growth  nehind  coneu  are  compared  with  ballistic 
range  data  taken  '-y  vlattary  (4)  and  •';lfcy  awl  -Sipwrly  (5)  and  Short  (Table  II). 

Ta  Section  2,  a  brief  resume  of  the  t.'.ec.i-»tical  relatior.'  is  include.i 
.a  order  to  provide  a  couvenieiii  refero-aee.  T-cv  .thock  shapes  for  t'v,  various 
todl=a  are  examined  in  Section  3*  It*  S'ictiv’'.  -  the  iflethods  used  to  calculate  the 
case  an-d  neck  properties  are  described  war  it  io.  .-•■i’C'iioo  p  tbs  coraplAtc  initial 
’rvipcii  vake  enthalpy  0i8triV.«<.ions  .are  given.  lascussion  of  the  aown°trecun 
t  urbulent  wake  prcpoi'tles  is  given  ie  Secr<,;-.ou  6  eau  conclusions  from  this  study 
in  Section  7* 

2.  SUMMARY  OF  T.TEORETICAL  AKAIY-Slg 

2.1.  Turbulent  Diffusion  of  Enthalpy  in  the  Inner  Wake 

Since  the  details  of  the  arAaiysis  used  here  have  been 
published  in  Ref.  (l),  only  essential  points  and  relations  are  listed 
below. 

1.  Howarth-Corcdnltsyn  variable  is  used  for  the  inviscid  regime 

where  y  =  y/d  aa.1  the  subs'cri-.jt.  L  iono^ies  4r.viBcid  or  outer  wako 

properties . 

2.  Inviscid  enthalpy  is  characterlred  by  e  two  jAi'jraia.  io.r  relation 
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g  (Y  )  «  Rg  (Y,) 


(8) 


m  =  0  for  two  dimensional  flow  and  m  =  1  for  axi-symmetric  flow’. 
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3.  Velocity  is  assumed  to  be  an  average  velocity  (u  =  ^  u^) 
where  0.8<p<0.9  so 
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4.  Howarth-Dorodnitsyn  variable  for  the  turbulent  core  is 
referenced  to  the  front  values  at 

5.  Profile  within  the  turbulent  core  is  characterized  by  two 
parameters 


(5) 


where  for  the  present  G 
following  will  be  used 


is  assiuoed  parabolic  and  the 
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6.  Reynolds'  hypothesis  of  similarity  between  the  turbulent 
transfer  of  mass^  momentum  and  energy  is  assumed.  Thus  only 
the  thermal  diffusivity  is  needed.  Townsend  (l)  found  that 
for  low  speed  wakes 

[uf  "  u(0)  ]l 


=  constant 


(6) 


where  is  a  measure  of  the  wake  width.  The  assumption  is 
now  made  that  the  hypersonic  wake  instantaneously  acts  like  a 
slice  of  low  speed  wake  or 


I 

where  K 


e  ^  =  1C  ^Uj  -  u(0)  jy^  ,  1  , 

can  be  related  to  Townsend's  IV,  by  X  = 

1  \-®f  *^1;, 


rl  (7) 


7.  Only  the  thcrin'3fiy«iamic  c-ficet  of  thr  s+s^tic  pressure  '/ari- 
ations  along  the  waKe  axis  is  conslaorea,  i.e.,  the  prosaurc 
graSient  terras  in  the  governing  equations  are  negiectert  tut 
the  equation  of  state  pA^ph  is  used.  By  iiitsg. -.ting  the  energy 
equation  across  the  wake  results  In  the  fo!.lowing  re.ictioa 
between  the  two  parameters  j^j  |'r| 


j^m+l 


(7-1) 
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Note  that  this  relation  gives  the  eemparison  of  ■inrt.uiear.eoas 
core  drag  to  initial  drag,  Thus  arymptoticslly  dowr  fiV,v*ati 
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Also  in  Eq.  8  the  fact  Is  used  that  the  turbulent  i.rx’i  Inrlscld 
Howarth  coordinates  (Y^,  Y^)  are  related  by 


~  Y„  Yy  ^  (constant)  (10) 
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8.  A  second  relation  between  3^^  (^)  and  is  obtained  by 

satisfying  the  energy  equation  along  tha  .vsiri.  llaaily,  using 
this  relation  and  Eqs.  (l  -  10)  leads  to  an  integral  equation 
for  the  wake  growth  Yj^  : 
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2.2  Turbulent  t'fass  Diffusion  in  the  Ituier  Wake 

0ne3  the  history  of  the  turbulent  diffusivity  along 

the  wake  axis  is  known,  turbulent  diffusion  of  a  particular  species 
in  the  inner  wake  can  bo  described  by  an  analysis  similar  to  that  for 
the  enthalpy.  It  is  assumed  the  mass  concentration  is  again  given  by 
a  two  rarameter  relation 


V'i'o 


Y  (~) 


} 


where  F  is  taken  to  be  gaussian.  The  result  as  given  in  Ref.  (l) 

for  the  variation  of  K  (5)  along  the  axis  is 
o  d 


(13) 


Once  the  mass  fraction  is  known  the  value  of  the  number  density 
per  unit  vol.ume  along  the  axis  is  calculated  from  the  relation 


The  relations  as  given  above  are  In  a  sll{^tly  different  form  from 
those  in  Ref.  (l).  This  is  because  the  over-all  invlBcid  drag  coefficient  Cjj 
was  introduced  as  a  normalizing  factor  in  the  former.  It  is  not  necessary  to*^ 
do  this  since  the  behavior  of  the  core  depends  only  on  the  location  of  the 
local  enthalpy  front  and  its  slope;  it  knows  nothing  of  the  over-all  drag  xmtll 
it  has  engulfed  the  complete  Inviscid  drag. 
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3.  BOW  SHOCK  SHAPE 


The  +,urhulent  wake  analysis  assiunes  that  the  outer  inviscid  enthalpy 
profile  is  invariant  with  axial  distance,  which  will  be  the  case  once  the  pressure 
returns  to  ambie;  .  Since  the  entropy  is  constant  along  any  streamline  in  the 
inviscid  flow,  only  the  shock  shapes,  including  the  wake  shock,  are  needed  to 
determine  the  desired  enthalpy  at  p  --  p^. 

The  general  configuration  of  the  inViscid  flow  field  as  calculated  is 
shown  in  Fig.  3  where  A  Is  a  constant  pressure  free  streamline  and  B  is  a  cylindri¬ 
cal  sting.  To  determine  the  wake  shock  actually  calls  for  knowledge  of  the  viscous 
base  and  neck  region,  since  these  determine  its  position  and  the  amoimt  of  recom¬ 
pression  necessary  for  the  flow  to  be  turned  back  parallel  with  the  axis  (angle 
of  A  and  position  of  B  in  Fig.  3)*  Strictly  speaking  the  bow  shock  is  also 
affected  by  the  base  region  through  the  shoulder  expansion  fan.  However,  the 
back  portion  of  the  expansion  fan,  which  is  sensitive  to  the  s;.oulder  expansion 
angle,  hits  the  bow  shock  far  downstream  where  the  shock  is  already  quite  weak. 

Thus  an  error  in  the  estimate  of  the  shoulder  expansion  angle  of  several  degrees 
was  found  to  make  negligible  differences  in  the  dowiibtream  bow  shock  shape,  and 
the  bow  shock  can  be  studied  separately  without  knowing  exact  details  of  the 
base  region.  The  latter  will  be  examined  before  the  second  shock  and  complete 
initial  wake  inviscid  enthalpy  profiles  are  given. 

For  each  of  the  configurations  studied  these  shook  shapes  were  obtained 
by  the  method  of  characteristics*.  Emphasis  is  placed  here  on  the  use  of  the 
shock  shape  rather  than  the  detailed  characteristics  solution  throughout  the 
flow  field  since  extension  of  characteristics  calculations  well  downstream  is 
apt  to  Introduce  significant  errors,  especially  in  high  entropy  gradient  flows 
(1,  6).  The  details  of  the  body  shapes  euid  flow  parameters  are  found  in  Fig.  8 
and  Tables  I  and  II. 

Blast  wave  theory  predicts  for  a  point  source  of  energy  that  the 
shock  should  be  parabolic  or 


*  In  the  ballistic  range  case  pictures  supplied  by  Dr,  Slattery  of  Lincoln  lab 
facilitated  in  checking  the  caiw-lrtlons. 
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Van  Hise  (7)  through  a  characteristics  study  demonstrated  that  for  a  hemis¬ 
phere  nose  at  high  ^ach  numbers  the  exponent  is  slightly  less  than  0.5.  Seifl’ 
(8)  and  V/hiting  point  out  from  experimental  range  firings  of  sphere-cones  that 
one  can  talk  about  three  shock  shape  regions,  two  of  which  can  be  approximated 
by  an  equation  of  the  form  of  Eq.  (l6),  but  each  with  different  exponents.  Con¬ 
sidering  even  more  detail  for  sphere-cones  if  the  body  is  sufficiently  slender, 
five  regions  cun  actually  be  said  to  exist.  This  is  seen  from  Fig.  4  in  which 
the  entropy  is  plotted  (or  shock  strength)  vs.  radial  shock  location  in  nose 
radii.  The  five  regions  are 

1)  high  entropy,  spherical  nose  segment  which  exists  until  the  last 
characteristic  from  the  sphere  intersects  the  shock 

2)  undershoot  and  recompression  around  initial  cone  section 

3)  conical  segment  if  the  body  is  long  enough 

4)  decay  starting  with  intersection  of  cone  shoulder  Mach  wave  and  shock 

5)  asymptote  downstream  where  shock  approaches  a  Mach  wave. 

In  Fig.  4  the  decay  regions  for  a  sphere  and  two  sphere  cones  l/lS 

are  shown.  The  relative  length  of  the  cone  segment  (or  nose/base  radius)  de¬ 
termines  when  the  expansion  fan  begins.  For  the  sphere  entropy  decays  smoothly 
from  the  nose  area.  For  Rjj/Rg  =  1/2  the  shoulder  expansion  fan  "catches"  the 
shock  at  the  limit  of  the  over  expansion  region  and  initiates  the  decay  there . 
For  Rjj/Rg  =  l/l2  the  shook  has  become  asymptotic  to  a  conical  shock  before  the 
7Ci  nose  radii  position.  If  the  body  were  an  ideal  pure  cone  the  shock  would  be 
constant  in  strength  until  the  shoulder  expansion  started.  The  axial  shock 
coordinate  corresponding  to  the  shock  position  is  written  along  the  curve. 

The  bow  shock  shapes  for  all  the  bodies  considered  have  been  plotted 
in  the  more  conventional  manner  of  _b  vs.  g  in  Figs.  ^  -  o.  Fig.  5  shows  the 
results  for  the  spherical  nose  bodies  (Cases  I  -  IV)  while  Fig.  6  is  for  the 
pure  cones  (Cases  A  -  D).  Kg.  7  includes  the  two  low  Mach  number  cones 
(Cases  E,  F)  as  well  as  that  for  =  22,  0  =  12°  (Case  D)  for  comparison. 
Finally  Kg.  8  compares  Case  II  with  a  ballistic  range  shot  by  Sommer  (9)  at 


Ames. 


In  Fig.  5  the  sphere  (Case  l),  which  is  the  same  as  used  hsjrore  in 
Ref.  (l),  gives  a  strai^t  line  with  slope  n^  =  .466^  which  is  close  to  the 
value  obtained  when  extrapolating  to  hi^.cr  Mach  numbers  from  the  data  of  Seiff 
(8)  and  Whiting,  The  sphere-cone  (Case  II)  has  a  break  in  the  curve  at  ^  =  3,7 
due  to  the  initial  effects  of  the  cone  section  .  Although  the  initial  segment 
(spherical)  of  the  sphere-cone  has  n^  =  .466,  the  latter  part  has  n^  =  .541. 

This  bears  out  a  fact  which  will  occur  consistently  here  and  was  mentioned  by 
Vagllo-Iaurln  (10)  and  Bloom,  namely  that  as  the  drag  coefficient  of  the  bodies 
becomes  smaller  this  decay  exponent  increases.  In  fact  for  Case  III  (sphere- 
cone,  Rjj/Rg  =  l/l2)  Uj,  =  .613. 

Case  III  serves  to  Illustrate  two  other  points.  First  is  the  far 

X 

downstream  approach  to  a  Msch  wave  which  begins  at  about  s  =  30*  Eventually  the 

o  ^ 

line  rfould  be  45  on  this  plot  and  n^  -  1.  This  region  is  also  seen  in  Case  IV 

which  is  the  compenrlson  of  a  sphere  fired  in  a  balllstlo  range  at  =  9»85 

Seiff  (8)  and  Whiting  and  present  characteristics  calculations. 

The  second  point  is  that  the  shock  becomes  nearly  conical  along  the 
body  and  remains  so  until  the  corner  expansion  hits  the  wave.  This  is  clearly 
seen  for  the  four  pure  cones  (Cases  A  -  D)  shown  in  Fig.  6. 

The  fact  mentioned  earlier  about  the  Increase  in  the  exponent  n  as  drag 
decreases  is  noticed  for  the  cones  here  with  «  .523  for  9  =  32°  increasing 
to  Hg  =  .613  for  9^  =  12”.  The  progression  of  this  slope  has  been  plotted  vs. 
cone  angle  for  the  cones  in  Fig.  9  and  vs.  drag  coefficient  in  Fig,  10.  It  is 
seen  to  be  neairly  linear  for  both  these  parameters.  Sphere  and  sphere-cone 
data  from  the  present  calculations  and  Refs.  (6)  and  (lO)  have  been  Included  in 
Fig.  10.  The  sphere  cone  exponents  for  bodies  with  drag  coefficients  between 
Cjj  =  0.27  and  0.48  seem  to  be  consistently  below  the  pure  cone  calculations, 
while  the  sphere  data  are  almost  in  line  with  the  cone  calculations. 


denotes  slope  in  nose  region,  n^  denotes  slope  in  wake  region  of  bow 
shock,  and  n^  denotes  wake  shock  slope. 

On  Figs .  5  6  the  short  vertical  line  along  the  curves  (marked  1 )  denotes 

the  shoulder  position  for  that  particular  body  while  the  dlagional  line  (s) 
denotes  position  of  first  shoulder  expansion  wave  intersection  with  the  shock. 
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Another  feature  which  will  turn  out  to  be  important  regarding  the  cone 

turbulent  wake  behavior  is  the  amount  of  mass  passing  through  the  constant 

strength  portion  of  the  shock,  i.e.,  the  position  of  the  shoulder  expansion  wave- 

bow  shock  intersection.  This  position  denoted  (Y.  has  also  been  plotted  in 

^se 

Figs.  9  and  10  and  it  is  seen  that  the  decay  exponent  Ug  and  behave  similarly. 

se 

At  very  large  angles  of  high  drag  coefficients  Y,  has  to  approach  an  asymptote 

^se 

because  Y,.  cannot  be  less  than  0.5.  The  point  at  =  I.5I  in  Fig.  10  is  from 
se 

a  60°  cone  calculation  taken  from  Vaglio-Iaurln  (10)  and  Bloom.  The  individual 
scattered  data  shown  in  Fig.  9  are  for  the  low-Mcich  number  experimental  data  luas. 

Fig.  7  shows  the  low  Mach  numoer  runs  for  the  pure  cones.  Again  n^ 
increases  as  cone  angle  decreases  and  also  as  the  Mach  number  decreases.  For 
M  =  5*9  and  9  =  12.5°  the  shock  Is  always  very  weedc,  in  fact  on  the  plot  it 
scarcely  deviates  from  a  45”llne.  This  case  gives  ng  =  .816. 

The  experimental  data  of  Sommer  (9)  shown  in  Fig.  8  was  the  only  experi¬ 
ment  available  on  a  slender  shape  at  high  Jtoch  nunibera.  This  body  was  fired  at 
*  17*95  and  is  compared  here  with  the  calculated  shape  for  the  sphere  cone  at 
=  22  =  1/2).  Although  the  experiment  was  at  a  slightly  lower  Mach 

number  and  9^  =  12.5°  it  is  felt  the  agreement  is  good. 

Finally,  all  the  shook  waves  computed  were  plotted  against  tlie  normal- 
ixing  factor  in  Fig.  11.  For  high  toch  nxaaber  the  downstream  decay  portions 
of  the  shocks  colfiipse  into  a  fan  all  passing  through  ^  =  1.26  and  -  2.0. 

In  the  nose  portion  of  the  shock  the  curves  deviate,  ifeicating  the  difference 
between  conical  type  growth  and  spherical  nose  shook  growth. 


+  For  slender  cones  the  Howarth  coordinate  in  the  Invlscid  wake  downstream  for 
a  given  enthalpy  level  is  very  nearly  equal  to  the  physical  coordinate  where 
that  particular  streamline  entered  the  bow  shock  since 

S  h 

Mass  =  p^u^n  y^g  =  2«d^  |  puydy  =  2j:p^d^  J  uYjdYj^ 

O  0 

and  at  p  ,  u  —  1.01  u  then  ^s  =  Y. . 

«  »  -T  I* 
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4.  BASE  AMD  MECK  HIOFERTIES 


Before  a  further  determination  of  the  inviscid  flow  field  it  :1s  neces¬ 
sary  to  examine  the  details  of  the  vJ.scous  hase  and  neck  regions.  This  was  doi*e 
only  approximately  since  no  complete  analysis  exists^.  It  was  assumed  the  flow 
is  as  postulated  hy  Chapman  as  shown  in  the  sketch  helow 


Sketch  1^  -  Base  and  Neck  Region 

with  a  recirculating  region  and  a  dividing  streamline  (A)  which  stagnates  on 
the  axis.  The  dividing  streamline  is  imbedded  in  the  viscous  portion  of  the 
flow  and  separates  the  recirculating  flow  from  that  which  passes  on  downstream. 
The  base  pressure,  free  shear  layer,  neck  size  and  initial  wake  pressure  are 
all  related.  Chapman  (ll),  Kuehn  and  Larson  give  the  pressure  rise  at  the  rear 
stagnation  point  in  the  viscous  portion  of  the  flow  along  the  dividing  stream¬ 
line.  This  pressure  rise  is  a  function  of  the  Inviscid  external  Mach  number 
(along  B)  and  the  velocity  w  =  j  along  the  dividing  streamline  (A)^^.  To 


Since  the  following  procedure  is  rather  detailed  and  some  of  it  is  accepted  to 
be  approximate,  the  q,ue6tion  arises  as  to  bow  much  does  the  downstream  wake  depend 
on  the  accuracy  of  the  Initial  quemtlties.  This  was  considered  somewhat  in  Ref. 
(l)  and  will  be  again  in  Section  6.4  for  tht  present  calculations. 

It  is  assumed  that  the  flow  is  brought  to  rest  isentropically  at  the  rear 

* 

stagnation  point  from  its  velocity  u  . 
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firvi  t'.ie  correct  base  nress  "c  tha  shoulder  expansion  angle  can  be  changed 

(vhich  then  changes  u  ana  thus  calls  for  an  iteration  procefiiare)  imtil  the  in- 

viscid  pressure  behind  the  shock  agrees  vith  the  viscous  value.  The  free  shear 

* 

layer  properties,  giving  u  ,  will  be  discussed  in  Section  4.1,1.  One  diffi¬ 
culty  in  carrying  out  the  above  procedure  in  axi-syirastric  flew  is  ubat  the 
inviscid  free  streailir^  is  not  straight  (see  Section  4.1.2),  thus  the  neces¬ 
sary  inviscid  pressure  ratio  to  turn  the  flow  parallel  to  the  wake  axis  varies 
depending  on  the  nacx  thicxness.  Although  the  free  shear  layer  properties 
were  examined  and  the  Chapsm  type  procedure  was  attenpted,  it  is  recognized 

that  only  a  rough  answer  could  be  obtained  because  of  the  approximations  neces- 
* 

sary  to  find  u  and  the  unknown  neck  thickiicss.  However,  once  acceptable  values 
for  these  quantities  are  found  along  with  the  pressure  behind  the  wake  shock 
properties  of  the  initial  wake  viscous  core  can  be  found  directly  (Section  4.2). 

4.1.  Free  Shear  layer 

We  now  consider  the  free  shear  layer  in  order  to  find  the 
* 

dividing  streamline  velocity  u  and  secondly  to  examine  the  effects 
of  curvature  of  the  inviscid  free  streainline  in  this  region. 

4.1.1.  Velocity  Along  Dividing  Streamline 

Chapman  (12)  has  analyzed  the  free  shear  layer 

problem  where  the  initial  boundary  layer  thickness  was  assumed 

to  be  zero.  For  this  case  he  finds  a  value  of  u  =  O.50. 

Denison  (13)  and  Bam  have  recently  re-examined  the  Chapman 

problem  with  an  initial  boundary  layer  thickness  and  find  cor- 
* 

rected  valuer  of  u  .  They,  iKwever,  do  not  take  into  account 
the  facts  that  the  body  has  a  finite  height  and  the  distortion 
in  the  free  shear  layer  profile  caused  by  the  expemsion  off  the 
body  shoulder.  Hamaitt  (l4)  did  the  same  type  of  analysis  in 
a  more  approximate  manner  assining  an  initial  cosine  shoulder 
profile.  His  final  results  which  are  given  in  Fig.  12  are 
used  here.  Fig.  12  gives  the  velocity  along  the  dividing  stream- 
line,  u  ,  as  a  function  of  the  Chapman  mass  parameter  Q  and  a 


r  = 


ux 


and  is  zero  for  the  present  along  the  dividing  streamline.  This 


can  be  considered  a  mass  injection  parameter  where  Qco  denotes  mass  injected 
resulting  in  u  <  u  (Q  =  0)  as  shown  la  Fig.  12. 


nixing  perarsster  g  vhere 


and  t  is  length  of  nixing  run,  6 .  .  is  the  initisl  trensfoiaed 

* 

incon^jressible  boundary  layer  thickness,  C  the  Chapeaan  Triscasity 

constant  ai^  the  Beynolds  nusber  is  the  local  free  shear  .layer 

value  based  on  1.  The  g  =  “  Unit  for  (g  =  O)  does  not  corres- 

pond  to  the  CSiapoan  value  of  u  ^  0.6  but  ratber  gives  s.  value 
* 

of  u^g  _  0.6o  due  to  the  approxLoations  made.  These  re¬ 

sults  can  then  be  used  to  find  a  correction  factor  to  be  applied 
to  the  Chapoan  value. 

Detailed  calculations  using  Eq..  (17)  and  Fig.  12  to 
determine  the  velocity  along  the  dividing  streamline  vere  made 
for  one  high  Ibch  number  case,  namely 

sphere  (Case  I,  Bef .  1)  u*  vas  taken  as  the  Chapman  value  =  0.6. 
As  mentioned  in  regard  to  the  vork  of  Denison  (13)  and  Baum  the 
initial  mixing  profile  after  the  shoulder  e^qpansion  off  the 
body  is  not  Blasius-llke  nor  a  cosine  type.  To  illustrate  the 
boimdary  layer  velocity  and  enthalpy  profiles  before  and  after 
the  shoulder,  as  -Talculated  for  Case  II,  are  shown  in  Fig.  13 
v'nere  the  profiles  have  been  e:qpanded  around  the  shoulder  by  an 
isentropic  stream  tube  method.  The  expansion  caiises  the  subsonic 
part  of  the  profile  to  conqjresS^with  the  result  of  the  flattened 
profiles  as  shown.  This  flattening  should  have  the  effect  of 
enhancing  the  Initial  mixing  and  hence  raising  u  .  For  the 
present  a  cosine  vas  approximated  to  the  expemslon  profile  of 
Fig.  12  with  the  result  that  u*  =  0.39  for  Case  II  at  }^  =  22 
and  ht  =  lOOK  feet  .•  To  find  the  Initial  profile  after  expansion 
the  a  priori  choice  of  expansion  angle  occurs  once  again,  but 
an  error  in  this  choice  of  several  degrees.. is  small  conqiared  to 
the  approximation  introduced  by  assuming  the  profile  to  be  cosine 
in  shape. 

Although  it  most  probably  is  not  correct  this  same  value  of  u*  was-assiuaed-- 
for  the  pure  cone  calculations.  The  effect  this  has  is  discussed  in  Section 
1».2.3. 
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4.1.2.  Inviscld  ?ree  Streaaiine 

ia  tvc  disensional  inviscid  flov  the  base  free  strean- 
line  esanating  frcn  the  shoulder  (A  in  Fig.  3)  is  strai^t,  tut 
for  axi-sycaetric  flow  this  streaaline  curves  inward  near  the 
axis. 


In  order  to  exaaine  this  point  the  shoulder  Frandtl- 
Meyer  ej^nsion  on  the  12°  pure  cone  was  set  successively  at 
initial  expansion  angles  of  0°,  -2°,  -  4°,  -  5°  with  the 
horizontal.  Fig.  14  shows  the  trajectories  of  these  stream¬ 
lines  compared  with  a  straight  line.  For  example,  for  9^  =  -  4° 
the  inclination  at  y/Bg  =  l/2  is  -  6°  and  approaches  -  10°  near 
the  axis.  These  angles  still  result  in  a  long  base  region,  but 
one  effect  of  interest  here  is  the  shock  strength  required  to 
deflect  the  streamline  back  parallel  to  the  axis.  As  shown  on 
the  figure  for  -  4°  the  difference  at  this  15ach  number 
(M^  =  17.3)  from  -  4°  to  -  6°  changes  p^^^  from  5  to  8. 

Thus  as  mentioned  before  the  size  of  the  neck  is 
now  important  in  determining  the  pressure  ratio  p^^^  . 

4.2.  Initial  Core  Properties 

From  Eqs.  (U)  and  (12)  it  is  seen  that  the  initial  viscous 
properties  of  the  core  which  are  needed  are  the  turbulent  front  lo¬ 
cation,  (Y^  )  ,  drag  in  the  core,  (C„  )  ,  the  axis  enthalpy, 

°fi  /i 


aid  the  initial  pressure  level. 


The  pressure  is  obtained 


from  the  characteristics  calculations  once  the  base  and  neck  details 
are  prescribed,  and  is  given  in  Section  5»1-  The  other  three  quanti¬ 
ties  are  related  through  Eqs.  (8)  and  (lO)  where 


1.^ .  X 

1^"  1  L  . 


®m+l  ” 


Remembering  that  p 


-6(constant). 


or  in  the  fora  of  Bq,.  (I2b).  Thus  it  is  really  necessary  to  calcu¬ 
late  independentlj'  two  of  the  input  quantities.  Two  procedures  were 

followed  depending  on  the  situation.  First,  where  a  picture  of  the 

+ 

neck  region  was  available  (Y  )  was  obtained  from  the  picture  and  inde- 

^f  i 

pendently  the  initial  drag  (C_  )  was  calculated.  Knowing  these  two 

®f  i 

resulted  in  a  value  of  h(0)^  from  Eq.  (l8)  or  (I2b).  However,  for 

most  cases,  in  the  absence  of  a  picture  (c„  )  and  h(0)  were  esti- 

“f  1  ^ 

mated,  and  a  check  was  made  by  Eq.  (18)  or  (I2b}  to  insure  that  a 

reasonable  value  of  the  physical  y^  occurred. 

■‘i 

4.2.1.  Initial  Core  Drag,  (C^  ) 

“f  i 

The  core  drag  is  assiuoed  to  be  made  of  two  parts: 
the  skin  friction  on  the  body  and  the  drag  associated  with  the 
pressure  rise  in  the  second  shook.  The  latter  is  found  easily 
once  the  wake  pressure  distribution  (Section  5.I)  is  found  (l) 
by  (P3  "  Pg)  axi-symmetric  bodies. 

Since  (C_  )  =  (C_  )  +  )  the  values  used  for  the 

“f  i  "f  skin  "f  1 

present  in  each  cate  for  (aC_. )  can  be  found  by  referring  to 

f  1 

Tables  I  and  II  where  all  values  of  (C_  )  and  (C..  )  are 

"f  i  ^f  skin 

listed.  Since  the  pressure  level  in  the  free  shear  layer  base 
region  for  a  blunt  body  is  usually  several  times  ambient,  where¬ 
as  the  slender  body  is  below  ambient,  )  for  the  blunt 

^f  i 

body  is  larger  than  for  the  slender  body.  For  example  the 
sphere  at  M  =  22  and  ht  =  lOOK  feet  has  {dC  )  =  .OO57  while 

*  4 


On  a  picture  the  physical  y_  ib  measured.  Section  4.2.2  discusses 
the  transformation  to  (Y,  ) . 

h  i 
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for  the  12°  cone  (aC^  )  =  .0018. 

i 

Skin  friction  for  pure  cones  as  a  function  of  cone 
angle  is  shown  in  Fig.  15^.  This  was  found  with  the  following 
relation  which  holds  for  conical  flow: 


i  °f'skin  ./T 


P®  u^ 


V  e,  R. 


(sin 


-1  cos  e 


1  -'B 


(’-9) 


h  2  ^ 

where  K  =  k"  .67 _ ^  .40  (Ref.  I5) 

Initially  there  were  two  questions  regarding  the 
cone  skin  friction.  First,  for  snail  cone  angles  the  cone  be¬ 
comes  very  long  and  (C^  )  should  become  quite  large.  The 

question  was  whether  this  occurred  at  cone  angles  of  practical 

interest.  Fig.  I5  shows  that  the  significant  increase  in 

(C_  )  begins  to  occur  at  very  low  angles,  in  fact,  below 
"f  skin 

one  degree  in  cone  angle.  Second  was  whether,  since  the  inviscid 
drag  is  much  less  for  a  slender  cone,  would  the  skin  friction 
now  become  larger  and  in  fact  dominant.  Considering  a  one  foot 
radius  body  and  angles  aroimd  10°  at  lOOK  feet  altitude  the 
laminar  friction  drag  is  still  small  compared  to  the  pressure 
drag.  If  the  body  were  small,  or  the  pressure  lower  (altitude 
higher)  this  would  no  longer  be  true.  For  a  body  radius  of  1/4" 
and  10°  cone  at  lOOK  feet  altitude  the  friction  is  over  one  half 
the  value  of  the  pressure  drag. 


4.2.2.  Turbulent  Front  Location 

The  relation  between  the  ^physical  location  and 
tran3fo:.*med  location  is  given  by  Eq.  (l2o) 


1  + 


(m+1)  H 


m+1 


i  I 


^f  '  m+1 

y  g(Y  )d  Y  L 

Lf  L^r 


(12c) 


The  pressure  drag  curves  shown  in  Fig.  IJ^will  be  discu^d  in  ^ction  ^.2. 

■’  ’  .  oc?  >;  » *  r* !. 
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I; 

so  if  the  physical  is  read  from  a  pictiire  f  is  found  directly 
by  iterating  to  get  a  solution  of  Eq,.  (I2c)  using  the  appropriate 
pressure . 

Pictures  of  the  neck  details  are  often  available  (as 
those  given  us  by  Dr.  Slattery)  for  ballistic  range  shots.  How¬ 
ever,  pictures  of  the  neck  details  for  higher  Mach  number  flows 
especially  for  the  slender  bodies,  do  not  exist.  Empirical 
correlations  of  neck  thickness  such  as  in  Ref.  (l6)  do  exist, 
but  in  general  they  are  not  reliable  enough  to  be  extrapolated. 
h.2.3.  Initial  Axis  Enthalpy  Level 

The  initial  value  of  the  axis  enthalpy  was  esti¬ 
mated  in  Ref.  (l)  by  assiuning  a  rear  stagnation  point  exists  at 
the  wake  neck  and  determining  the  stagnation  enthalpy  on  the 
stagnated  streamline.  This  stagnated  streamline  is  the  dividing 
streamline  considered  previously.  The  recompresslon  of  this 
streamline  to  zero  velocity  on  the  axis  is  assumed  to  occur 

rapidly  without  any  energy  loss.  To  find  the  stagnation  energy 

* 

along  the  dividing  streamline  it  is  necessary  only  to  find  u 
since 

^  s  *  \ 

1  -  _S_  u  4  _±_  (20) 

^  I 

where  H  is  local  stagnation  enthalpy  and  h  is  the  known  wall 

* 

enthalpy.  For  a  cold  wall  Chapman's  value  of  u  would  give 
h(0)j^  =  At  hypersonic  speeds  this  is  indeed  quite  hot. 

This  value  was  taken  for  the  sphere  as  calculated  in  Ref.  (l). 

For  the  case  calculated  in  detail  here.  Case  II,  in  which 

u*  =  0.39  (Section  4.1.1),  \  was  taken  as  0.13  thus—  =  h(0)^/H^ 

n  H 

“ 

=  0.47  from  Eq.  (20). 

Rather  than  repeat  thi.s  involved  procedure  for  the 
pure  cones  the  initial  h(0)^/h^  value  was  arbitrarily  taken  to  be 
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thio  caii'.o  value,  l.o,, 

o  bly  not  correct,  but  an  shown  later  it  turns  out  that  the  cor.e 
equilibrium  wake  behavior  is  relatively  insensitive  to  the  in¬ 
itial  enthalpy  level.  But  for  the  non-equilibrium  case  an  exact 
dctcrmlnutlon  of  thin  enthalpy  level  for  a  slender  body  is  still 
a  moat  imixjrtant  question. 

For  the  experimental  range  firings  it  v.’as  not  neces¬ 
sary  to  determine  h(0)j^  independently,  since  usually  the  direct 
method  of  measuring  the  nock  thickness  from  a  picture  was  used 
and  h(0)^  found  from  Eq.  (l8).  However,  for  these  cases  values 
of  h(0)j,/1l^0.5  always  occurred. 

4.3.  Wake  Shock 

With  the  shoulder  expansion  angle  and  the  neck  thickness 
determined  the  wake  shock  can  now  be  inserted  in  the  Inviscid  flow. 

To  review,  the  characteristics  calculations  were  made  for  the  geometry 
shown  in  Fig,  3  where  A  is  a  constant  pressure  free  streamline  and  B 
is  an  inserted  cylindrical  sting.  The  resulting  wake  shock  shapes 
are  shown  in  Figs.  5  and  7*  These  include  both  the  high  Mach  number 
results  and  those  for  the  experimentally  fired  cones.  The  nose  is 
so  small  for  =*  that  the  woke  shock  is  taken  to  be  the  same 

as  the  pure  cone  (Case  D),  Initially  each  of  these  wake  shocks  travel 
outweird  quite  rapidly  except  for  the  latter  cases  (III  or  D)  which 
seem  to  take  some  time  for  the  shock  to  strengthen.  The  wake  shock 
approaches  the  position  of  the  bow  shock  with  eventually 

becomes  aLnost  parallel  with  it  as  both  decay  to  a  Mach  wave  together. 

5.  PRESSURE  DISTRIBUTION,  PRESSURE  DRAG  AMD  INITIAL  WAKE  EWTHALPY  DiarRIBlglONS 
,^.1.  Pressure  Distribution 

The  pressxire  distributions  in  the  nose  region  and  along  the 
afterbody  for  each  of  the  bodies  calculated  are  shown  in  Figs.  I6-I8. 
Fig.  16  shows  the  four  sphere,  sphere-cone  cases  at  =  10  and  22, 

Fig.  17  the  pure  cones  at  =  22  and  Fig.  I8  the  two  experimental 


=  46  at  =  22.  This  is  most  proba- 
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cone  cases.  The  axial  coordinate,  used  in  the  figures  was  measured 
from  the  apex  of  the  shock  in  order  to  keep  the  coordinates  consistent 
with  those  for  the  shock  shape.  For  che  most  part  both  a  cylindrical 
and  wake  afterbody  are  shown.  Again  Refs.  7>  10  and  1?  present  detailed 
pressure  distribution  studies  for  air  so  only  a  few  comments  as  to  the 
present  particular  cases  will  be  made  here.  First,  if  the  flow  at  the 
body  shoulder  is  expanded  to  be  parallel  to  the  free  stream^  the  re¬ 
sulting  pressinre  is  not  much  different  from  that  which  will  occur  be¬ 
hind  the  wake  shock  (where  the  flow  is  again  parallel  to  the  free 
stream).  This  is  expected  since  the  wake  shock  is  relatively  weak. 

Thus  at  M  =22  the  initial  wake  pressure  for  the  blun+-  body  is  S'  24 

while  for  the  slender  cone  ^  ^  1  -  2.  This  is  one  of  the  fundamental 

Pa, 

differences  between  the  hypersonic  blunt  body  and  slender  body.  In 
general,  the  pressure  behind  the  wake  shock  for  the  blunt  body  is 
quite  high  whereas  for  the  slender  body  it  is  nearly  ambient. 

Secondly,  the  pressure  with  the  wake  shock  included  rapidly 
approaches  the  distribution  which  would  exist  for  the  same  forobody 
but  with  a  cylindrical  afterbody  (instead  of  wake).  In  fact  by  10 
diameters  behind  the  wake  shock  the  pressure  is  very  near  to  the 
cylindrical  value 

As  Illustrated  in  Fig.  17  and  pointed  out  by  Vaglio-Iaurin  (lO) 
and  Bloom,  the  cone-faced  bodies  at  high  Mach  number  exhibit  a  down¬ 
stream  recompression  on  the  cylindrical  afterbodies,  which  exists  con¬ 
sistently  for  all  cone  angles  considered  here.  This  also  is  seen  from 
Fig.  17  to  exist  when  the  wake  shock  is  inserted  (17) • 

The  calculations  for  the  two  low  speed  experimental  runs  eire 
shown  on  Fig.  18. 


*  For  a  sphere  it  already  is  parallel  to  the  free  stream  at  the  shoulder. 

This  is  not  true  regarding  the  temperature  since  any  entropy  rise  persists 
in  an  Inviscid  flow. 
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5 • 2 .  Invlscld  Forebody  Drag  Coefficient  (Pressure  Drag) 

The  inviscid  drag  coefficient  for  pure  cones  as  a  function 
of  the  cone  angle  is  shown  in  Fig,  15  and  for  various  sphere-cone  coiiibi- 
nations  (from  Ref.  (I8))  in  Fig.  19* 


Previously  (l)  the  drag  coefficient  for  a  sphere  had  been 

taken  as  unity.  The  more  detailed  characteristics  results  indicate 

here  a  value  of  =  O.94.  Characteristics  calculations  made  present- 
Q  o 

ly  for  the  12  cone  are  shown  as  the  circled  line  in  Fig.  19»  The 
difference  for  this  one  angle  is  that  the  calculations  were  carried 
very  far  downstream  such  that  a  pure  cone  was  approached  (Rjj/Rg  =  *075)' 

The  approach  to  the  slender  cone  shows  that  the  drag  coefficient 
actually  undershoots  the  cone  value  by  a  small  amount  (19)>  and  that 
for  a  value  of  as  high  as  0.2  the  drag  coefficient  is  almost  the 

same  as  for  a  pure  cone.  This  is  due  to  the  pressure  undershoot  on 
the  cone  surface  when  a  small  spherical  nose  is  put  on  the  body. 

Thus  it  is  necessary  to  blunt  the  nose  a  surprising  amount  before 
any  significant  effect  on  the  over-all  drag  coefficient  occurs.  This 
fact  will  have  significance  later  in  determining  the  wste  properties. 


5. 3*  Initial  Wake  Inviscid  Enthalpy  Distribution 

The  inviscid  enthalpy  distributions  downstream  at  p  =  p- 

00 

are  shown  in  Fig.  20  for  M  =  5»9  and  Figs.  21  and  22  for  M  =  22. 

00  00 

The  profiles  as  shown  are  in  the  liowarth  coordinate  which  is  used  in 
the  analysis. 


For  the  case  of  M  =5.9  the  enthalpy  both  with  and  without 

40 

the  wake  shock  is  shown  in  Fig.  20.  The  bow  shock  profile  is  constant 

until  the  expeuislon  fan  occurs.  The  sudden  change  in  slope  here  plays 

an  important  role  in  the  wake  history.  The  wake  shock  contribution 

is  significant  near  the  axis  where  the  flow  deflection  is  the  greatest. 

However,  as  Indicated  by  the  line  labeled  (Y^  )  most  of  this  peak  is 

^f  i 

engulfed  in  the  viscous  core  at  the  neck.  The  effect  of  the  second 


shock  is  to  "smear"  out  the  vertical  front  presented  by  the  bow  shock. 

For  this  low  Mach  number  even  the  peak  inviscid  enthalpy  at  =  0  ?s 

h  “ 

relatively  "cold",  giving  only  =  1.2.  This  effect  tends  toward 

CO 

computational  difficulties  when  computing  quantities  such  as  drag 
which  depend  on  very  small  differences. 

The  high  Mach  number  cone  shown  in  Fig.  21  ir.dicates  again 
that  most  of  the  "hot"  inviscid  region  due  to  the  wake  shock  should 
be  already  swallowed  by  the  "neck".  In  this  case,  to  the  accuracy  of 
the  computations,  the  wake  shock  gave  no  noticeable  contribution  be¬ 
yond  the  knee  in  the  curve  (Y,  >Y,  ).  For  this  cone  even  at  M  =22 

and  including  the  wake  shock  the  maximum  inviscid  enthalpy  (g  ^  3)  is 

still  "cold"  compared  to  the  viscous  core  value  used  here  =  k6), 

h 

00 

Primarily  shown  on  Fig.  22  is  a  compilation  of  enthalpy  pro¬ 
files  for  M  =  22  Including  both  families  of  sphere-cones  and  pure 
cones.  For  Rjj/Rg  =  i/lS  the  extremely  small  hot  pulse  (normal  shook 

enthalpy  expanded  to  p  =  p  )  occurs  because  of  the  small  nose  radius, 

00 

The  overshoot  on  the  cone  region  of  the  shock  shape  as  discussed 
earlier  is  evident  here  by  the  dip  in  the  enthalpy  profile,  followed 
by  the  approach  to  the  cone  value  at  Y.  =  0.*f,  Boundary  layer  calcu- 
lations  at  the  cone  shoulder  show  that  the  mass  in  this  inviscid  hot 
pulse  is  entirely  engulfed  in  the  boundary  Is.yer  by  that  point,  thus 
it  is  Incorrect  to  carry  it  into  the  wake  calculations.  In  fact,  as 
shown  on  Fig.  22  the  enthalpy  undershoot  is  also  within  the  viscous 
region  at  the  time  the  flow  reaches  the  neck.  Thus  if  the  nose/base 
radius  ratio  is  as  small  as  l/l2  pure  cone  flow  can  be  used  for  equi¬ 
librium  wake  calculations. 

For  Rjj/Rg  =  the  width  of  the  inviscid  enthalpy  pulse 
has  been  reduced  from  the  sphere  by  a  factor  of  2-3  which  is  certainly 
predictable  from  the  drag  coefficients.  Here  the  viscous  front  is 
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still  within  a  very  hot  region  but  near  the  point  of  rapid  change, 
'i'his  will  effect  the  nature  of  the  wake  history  as  seen  later. 

Besides  the  12°  cone  Fig.  22  shows  pure  cone  enthalpy  pro¬ 
files  for  6^  =  22°  and  32°.  From  these  curves  comparing  the  general 
nature  of  the  sphere  cones  and  pure  cones  it  can  be  seen  how  differ¬ 
ently  the  nose  shapes  distribute  the  waKe  enthalpy.  One  very  vivid 
comparison  is  the  22°  cone  and  sphere-cone  (S^  =  12°, 
which  have  the  same  over-all  drag  but  from  Fig.  22  are  seen  to  be 
very  much  different  in  distribution. 


In  conclusion  the  inviscid  flow  field  can  be  affected 
drastically  by  shaping  the  body.  We  now  turn  our  attention  to  the 
calculation  of  the  downstream  wake  properties. 


6.  DISCUSSION  OF  DOWHSTREAM  WAKE  PROPERTIES 


In  order  to  find  the  downstream  wake  properties  Eq.  (U)  along  with 
its  component  parts  lias  been  programmed  on  a  computer.  The  necessary  inputs 


arc  the  initial  values  of  (Y,  )  ,  (C„  )  , 

^f  i  °f  i 


li 


along  with  tho  downstream 


pressure  history,  2  vs.  and  the  radial  enthalpy  distribution  g(Y- )  vs.  Y^ . 


The  enthalpy  distribution  was  put  in  numerically  as  determined  directly  from 
the  shock  shape.  Previously  (l)  the  distribution  was  approximated  analytically 
in  such  a  manner  as  to  make  the  total  drag  correct.  The  advantage  of  using  the 
numerical  value  is  that  no  distortion  of  the  profile  occurs  in  the  sensitive 
inner  portion  by  forcing  an  over-all  curve  fit.  Also,  if  only  partial  shock 
shapes  are  known  the  wake  can  be  found  at  least  for  some  distance  downstream. 


A  question  occurs  concerning  the  initial  eixial  station  at  which  mixing 
calculations  should  be  started.  This  occurs  because  the  analysis  does  not  in¬ 
clude  a  pressure  gradient  and  the  pressure  drops  initially  very  rapidly,  es¬ 
pecially  for  a  blunt  body.  The  procedure  adopted  in  Ref.  (l)  was  again  used, 
namely  it  was  arbitrarily  assumed  no  mixing  occurs  until  the  pressure  is  down 

to  o/  =4.  The  wake  is  assumed  to  simply  expand  by  the  square  root  of  the 
*  P 

pressure  ratio  from  the  neck  to  the  Initial  mixing  station,  i.e.,  v/here  p/  =  4. 
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Since  for  a  blunt  boay  at  high  speeds  the  pressure  behind  the  shock  is  about 
24  p^  this  results,  in  practice,  in  starting  the  mixing  calculations  6-8  di¬ 
ameters  downstream.  For  a  slender  body,  however,  (6^ps12°)  the  calculations 
are  started  right  at  the  neck  since  the  pressure  is  less  than  4  initially, 

6.1.  Wake  Growth 

The  three  features  brou^t  out  below  deal  with,  1.)  shape 
variation  comparing  wakes  of  blunt  and  sharp  bodies,  2.)  shape  vari¬ 
ation  effects  on  bodies  with  the  same  drag  coefficient,  and  3«)  alti¬ 
tude  (Reynolds  number)  effects, 

6.1.1.  Effect  of  Shape  Variation  (Blunt  vs.  Sharp) 

The  available  experimental  wake  growth  data  are 
from  Slattery  (4)  and  Clay  and  Siperly  (5)  and  Short.  A  com¬ 
parison  of  the  theoretical  and  measured  growth  for  the  bal¬ 
listic  range  data  of  Slattery  (4)  and  Clay  is  shown  in  Fig.  23. 
The  figure  also  shows  their  measurements  for  spheres  and  the 
coirespondlng  calculations  made  in  Ref.  (1).  The  difference 
in  character  in  the  sphere  growth  and  cone  growth  can  be  seen. 
The  dotted  line  is  the  asymptotic  value  for  the  cone  calcu¬ 
lations  and  it  falls  below  the  sphere  calculations  by  the  cube 
root  of  the  drag  ratio.  The  present  theoretical  line  for  the 
cones  is  on  the  low  side  of  the  data  and  is  lower  than  the 
scatter  from  100  <  3/1,000,  The  theoretical  inflection  point 

X  ^ 

around  ^  =  600  is  the  point  at  which  the  front  breaks  out  into 
flow  coming  through  the  expansion  region  of  the  shock  and  the 
drag  (as  small  as  it  is)  is  engulfed  rapidly.  It  would  appear 
something  like  this  happens  around  ^  =  100  in  the  experiments. 

The  comparison  with  the  Siperly  (5)  and  Short 
data  is  shown  in  Fig.  24.  The  data  was  taken  at  NASA-Ames. 

This  body  is  a  blunter  cone  and  at  slightly  higher  velocity 
than  the  Lincoln  data.  The  data  indicates  a  rather  sudden  in¬ 
crease  in  wake  width  around  3  r,100.  This  is  the  typical  be- 

d 

havlor  predicted  for  a  cone;  however,  the  "explosion"  predicted 
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from  the  theory  occurs  later  at  around  ^s»300.  The  early  humps 
in  the  curves  around  ■:  10  occur  due  to  the  nressure  decrease 

X  ^ 

to  ambient.  At  -a:l,000  most  of  the  drag  should  be  engulfed  by 
a 

the  core  and  the  curves  should  approach  their  asymptotic  behavior. 

Growth  curves  for  three  cases  at  high  :iich  number 

arc  shown  in  Fig.  25  to  illustrate  further  the  conti'ast  between 

o 

blunt  and  si.cndor  bodies.  Those  shown  here  are  the  12  cone, 
sphere-cone,  sphere  fami3.y  (Cases  I,  II  and  D).‘ 

Here  again  the  cone  curve  is  significantly  below 
the  other  shapes.  The  sudden  explosion  at  about  j*“30  for 
Rjj/Rjj  =  l/2  is  due  to  the  previously  mentioned  fact  that  the 
front  started  in  a  hot  region  but  very  close  to  a  rapid  drop 
off  in  the  enthalpy.  The  cone  wake  growth  is  near  a  cube  root 
throughout  except  around  jsaTCO  where  it  hits  the  sudden  drop 
in  inviscid  enthalpy  as  seen  in  Fig.  21.  The  humps  in  the 
growth  curves  for  low  values  of  j  (<  50)  are  again  due  to  the 
pressure  decay  characteristics. 

The  rate  of  turbulent  core  growth  depends  on  the 
amount  of  Inviscid  drag  swallowed.  In  fact,  over  the  whole 
regime  a  rough  qualitative  estimate  (l)  states  that  for  an 
axi- symmetric  wake 


(21) 


If  Eq.  (21)  is  nonnalized  by  ^  C^A,  where  A  is  the  base  area, 
then  (2) 


The  only  difference  from  a  low  speed  wake  is  that  C  (y  )  i& 

Df  T^. 

a  changing  function  except  in  the  constant  conical  region  for 
a  cone  and  the  far  downstream  asymptote  where  (Y  )  /C  =1.0 
(See  Eq.  $). 

For  the  early  stages  of  the  cone  growth,  rather 
than  an  approxinate  equation  of  the  form  of  Eq.  (22),  the  full 
Eq.  (ll)  can  be  simplified  considerably.  As  long  as  the  turbu¬ 
lent  front  for  the  cone  is  growing  into  a  uniform  region,  which 
is  "cold",  no  relative  drag  compared  to  the  initial  front  is 
being  engulfed.  For  the  cone  up  to  the  point  where  the  front 

u 

begins  to  decrease  g(Y^ )  =  1,  f  =  1  +  H  and  from  Eq.  (l2o) 

^  h" 


1  +  H 


Also,  since  the  slope  of  the  front  is  zero  and  no  relative  drag 
has  been  engulfed,  =  Fg  =  1.  With  these  facts  Eq»  (ll)  be- 


The  axis  enthalpy  is  now  =  (l+K)  +  >  so  Eq.  (ll') 

t'-l  6 


becomes 


J  (Yj^“'"^  6  )  (l+H)  + 


const . 

Y  6 
^f 


1 

m+1 

=  const: 


\  m+1 

d^i“ 


The  axis  enthalpy  term  in  the  ^  j  brackets  has  the  effect  of 
slowing  down  the  initial  growth,  but  becomes  small,  fairly  soon 


and  when  Y^^  6,  Y, 


(C_  )  •  ^  similar  to  Eq.  (21). 

Uf  i  a 


Thus  as  h(0) _ ^  h^  and  prior  to  reaching  the  break  in  the  enthalpy 

curve  the  cone  growth  should  go  as  the  cube  root  and  be  nearly  pro¬ 
portional  to  (C^  )  .  This  is  shown  for  Y^  in  the  sketch  below 
®f  i  % 


-^hlcl.  is  ti- 
The  physical 
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X 
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calcaLitioss  for  the  12°  cooe  at  Zi  =  22. 

A 

shown.  S:e  ^•aryin2  difference  cetveea 
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Sketch  2.  Initial  Cone  Vake  Growth 

At  the  far  downstreaa  asynptote  where  all  the  drag 
has  been  engulfed  by  the  turbulent  core  Eq,.  (il)  becoaes 


y  ^ 


K  G  (0) 


4=‘^^G 


(C  )  ^ 
'  d'  d 


(23) 


n+1 


which  when  norraalited  with  ^/CjjA  becoaes  like  Eq.  22  where 

C„ 


{X,  )/c„  =  1,0 
3 


“Df  D 


Const,  —z. 


(24) 


Thus  plotting  the  wake  growth  nornalised  by  ^^C^A  gives  one 
curve  asymptotically  for  all  bodies  regsjrdless  of  shape  and 
size.  This  is  shown  in  Fig.  26. 

6.1.2.  Effect  of  Shape  Variation  (Same  Drag  Coefficient) 

The  blunted  cone  Rj/Rg  =1/2  and  the  22°  pure  cone 
have  the  sajDe  drag  coefficient  {C^  =  .30).  The  inviscid 
enthalpy  profiles  were  seen  in  Fig°  22  to  have  campletely  differ¬ 
ent  enthalpy  distributions  although  the  drag  Integral  is  the  same. 
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^rcs  the  discussion  in  tbe  previous  section  on  general  blunt  vs, 
slender  vale  growth  one  dgjt  expect  tbs  early  portions  of  tbe 
wake  of  these  two  bodies  to  be  different  and  this  is  seen  to  be 
true  in  rig.  2?.  Eventually  the  curves  coincide  but  until 
gajljOCO  tbe  behavior  for  these  two  todies  of  tbe  sane  drag  (ana 
area)  is  cuite  different. 

6.1-3-  Effect  of  Altitude  Variation 

The  effect  due  to  altitude  is  caused  prlaarily  by 

tbe  Heynolds  rubber  tbrou^  (C^  )  .  In  coeparing  snail  scale 

f  i 

ballistic  range  data  with  large  body  hi^  Ifecb  conditions  this 

effect  can  be  very  irnortant.  To  illustrate,  three  hi^  tbcb 

cases  were  calculated  at  altitudes  of  60K,  lOOK  and  IpOK  feet. 

For  the  higher  altitudes  or  lower  Reynolds  number,  the  skin 

friction  becomes  a  larger  percentage  of  the  total  drag.  This 

results  in  the  turbulent  core  starting  farther  out  on  the  in- 

viscld  front,  which  In  turn  tends  to  reduce  the  Initial  growth 

dlffersaces  seen  to  exist  previously.  Indeed,  if  (C^  )  is 

^f  i 

large  enou^,  all  that  occurs  is  tbe  asymptotic  cube  root 
growth. 

The  importance  of  this  effect  is  seen  in  Fig.  27, 

where  curves  for  both  lOOK  and  IpOK  feet  are  included,  and 

Fig.  28  for  the  12°  cone.  In  Fig.  27  the  difference  between 

the  two  bodies  is  smaller  at  I50K  feet  than  at  lOOK  feet.  In 

fact,  in  Fig.  28  for  the  12°  cone  at  I50K  feet  no  hrap  in  the 

curve  is  observed  because  the  front  starts  outside  of  the 

enthalpy  break.  The  reason  that  the  6OK  feet  curve  is  above 

that  for  lOOK  feet  is  because  turbulent  skin  friction  was 

assumed  to  exist  at  this  altitude.  If  transition  on  the  body 

occurs  at  an  altitude  different  from  what  was  assumed  here, 

obviously  (Cn  )  changes  and  consequently  also  the  downstream 
i 

wake. 
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StlUstlc  riiK-c  art:  sade  with  fairly 

ii-  iiii  hc-iaa  ct  ratucci.  pressure  so  cars  laist  ce  cxcrcisea  in 


•y  SiSase  var_- 


caoosin^  the  ca.xirii-jcnts  if  one  is  iooi';in5  for 


6.2.  Snthal:^  and  Valogity  Sistribations 

Perhaps  the  aost  cienificant  effect  of  shapinc  can  bo  seen 
fros  thv.  cntlni^  anC  velocity  distrihutiors;  as  shown  in  Tigs.  29 
and  30-  Fig.  29  shows  the  enthalpy  and  velocity  distributions  aioi!g 
the  axial  strcanline  and  the  front  for  the  12°  cone  -  sphere  fasaly 
(Cases  J,  II  and  D).  These  arc  plotted  against  the  "reduced"  co¬ 
ordinate  ^  — -  Since  the  analysis  assures  a  constant  stagtULtion 
v%*‘ _ 

enthalpy  u  ■^j  H  -  h,  and  the  velocity  deficit  along  the  axial  streac- 

-  u(o) 

line  noraalized  by  the  free  streaa,  -  ,  is  shovm  •  n  the  rij^.t 

CO 

hand  ordinate.  Fig.  30  shows  the  enthalpy  and  veloeity  distributions 
along  the  axial  streaaline  at  two  altitudes  for  the  two  bodies  with 
the  sane  drag  (Cases  II  and  B). 

It  is  noted  that  reducing  the  nose  radius  significantly  re¬ 
duces  the  wake  temperatures,  especially  for  the  cone  which  drops  an 
order  of  aagnitiide  faster  than  the  sphere.  In  fact  the  cone  drops 
so  rapidly  it  appears  the  very  early  history  is  insensitive  to  the 
exact  initial  enthalpy  level  chosen.  This  general  behavior  is  not 
evident  beforehand  because  the  front  spreading  rate  is  much  slower 
for  the  cone.  Conceivably  the  low  cone  wake  growth  rate  could  cause 
any  hi^  core  temperatures  to  persist  far  downstream.  The  long 
persistence  cf  a  high  core  temperature  is  enhanced  by  and,  in  fact,  de¬ 
pends  on  the  core  engulfing  hot  fluid  or  drag  from  the  outer  region. 

For  an  axi-syicnetric  body  the  enthalpy  excess  of  the  axis  above  the 
front  is  given  by  Eq.  (8)  or  (l8) 

h(0)  -  h^ 
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where  C  is  the  Initial  plus  engulfed  drag  in  the  core.  If  the  lao- 

nentun  defect  relative  to  the  front  increases,  then  Cjj  increases  with 
Y-  thus  keeping  (h(O)-h  )  large.  For  the  sphere  this^occurs  such 

that  (h(O)-h^)  is  aljnost  constant  for  hundreds  of  diasneters,  and  the 
overall  temperature  level  drops  only  as  the  front  temperature  de- 
creeises  which  is  relatively  slow.  For  the  cone  practically  no  invlscid 
drag  relative  to  the  initial  front  is  engulfed  for  hundreds  of  diameters 
(the  inviscid  enthalpy  profile  is  almost  constant),  thus  any  core  ex¬ 
pansion  at  all  has  directly  the  effect  of  drastically  reducing  the 
temperat'orc- 

The  velocity  follows  a  similar  history.  For  the  cone  the 
velocity  is  onl’r  5/^  different  from  free  stream  by  125  units  of 

- —  (—30  boily  diameters),  whereas  for  the  sphere  it  takes  nearly 

650  units  of  — - (—600  body  diameters).  Another  interesting  com- 

parison  deals  vith  the  velocity  slowdown  at  some  level  of  radar 
frequency,  say  UHF  (n^  (critical )a.1.6*10^  electrons/cm^).  At  this 
point  for  the  sphere  a  change  in  velocity  of  100  ft. /sec.  would  occur 
in  a  distance  of  dix.r~2^  feet  (assuming  d  =  2  feet)  while  for  the  cone 
Ax~l  foot. 

In  comparing  the  enthalpy  and  velocity  histories  of  the 

two  bodies  witii  the  same  drag  (Fig,  30,  Cases  II  and  B),  it  is  seen 

they  are  quite  different,  as  were  the  woke  growths.  The  22°  cone 

o 

follows  more  closely'  uhe  drop  of  the  12  cone  rather  than  the  blimt 

body.  Also,  as  in  the  wake  growth,  the  effect  of  increasing  (C„  ) 

"f  i 

(raising  the  altitude)  decreases  the  difference  between  the  bodies. 

6.3.  Electron  Decay 

6.3,1.  Equilibrium  Flow  Conditions 

The  effect  of  the  drop  in  enthalpy  with  decrease  of 
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nose  bluntness  is  felt  directly  on  the  electron  decay  as  showi 
in  Fig.  31*  This  shows  the  axial  electron  density  for  the 
cone,  sphcre-conc  and  sphere  (Cases  D,  I  and  II).  Both  the 
cone  and  sphere-cone  electron  concentrations  drop  rapidly  in 
comparison  to  the  sphere. 

The  depen:\ence  of  the  equilibrium  electron  concen¬ 
tration  on  drag  coefficient  for  this  class  of  bodjes,  i.e-, 
o 

sphere-12  cone  series,  can  be  seen  by  cross  plotting  the 

distance  it  takes  to  drop  to  various  electron  levels  vs,  the 

drag  coefficient.  This  is  shown  in  Fig.  32  where  the  electron 

levels  have  been  arbitrarily  chosen  as  those  corresponding  to 

critical  frequencies  for  C  band  (56OO  me),  L  band  (13OO  me), 

and  UHF  (400  me).  The  surprising  result  occurs  that  the  decay 

length  is  linear  with  the  drag  coefficient.  In  view  of  the 

other  results  concerning  bodies  with  different  cone  angles 

this  curve  cannot  be  applied  to  any  general  shape.  It  should, 

o 

however,  hold  for  other  sphere-cone  bodies  with  a  12  cone 
angle. 

6.3.2.  Rire  Mass  Diffusion 

The  electron  decay  for  pure  diffusion  is  also  shown 
in  Fig.  31  for  the  cone  and  sphere.  This  calculation  makes  use 
of  Eqs,  (14)  and  (15)  and  assumes  a  foreign  species  (e.g.,  ab¬ 
lation  products)  to  be  deposited  into  the  wake  through  the  neck 
with  an  initial  gaussion  radial  distribution  and  subsequent 
zero  recombination.  The  diffusivity  used  is  that  which  is 
found  for  the  enthalpy.  Fig.  31  shows  this  type  of  diffusion 
for  the  cone  to  be  even  slower  than  for  the  sphere .  Since  no 
recombination  takes  place,  the  slow  diffusion  results  from  the 
slow  cone  wake  growth.  Very  far  downstream,  after  the  "explosion" 
in  cone  growth  the  rate  will  eventually  be  the  same. 

Except  for  the  pure  diffusion  calculations  of  this 
section  the  assumption  of  thermodynamic  equilibrium  has  been 
made  everywhere  in  this  report.  This  means  that  as  far  as  an 
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application  to  the  re-entry  prohlem  is  concerned  the  results 
will  not  hold  above  a  certain  altitude,  which  is  approxJ.mately 
lOOK  feet  (see  Ref.  2)  due  to  non-equilihrium  effects.  The 
possible  importance  of  non-equilibrium  effects  is  emphasized 
by  the  results  cf  Fig.  31  in  noting  the  difference  between  the 
equilibrium  curve  on  the  one  hand  and  the  ultra-conservative 
zero  recombination  pure  diffusion  on  the  other  band. 

6.4.  Effect  of  Accuracy  of  Initial  Wake  Conditions 

An  effect  of  inaccurately  choosing  the  initia.,  conditions 
for  the  wake  can  be  inferred  either  from  the  enthal.py  histories  or 
the  wake  growth  curves  in  Fig.  27-  For  the  sphcrc-ccne  =  l/2) 

the  difference  in  (C^  )  at  lOOK  feet  and  I5OK  feet  is  almost  exactly 
*^f  i 

a  factor  of  two.  Beyond  the  relative  differences  in  wake  behavior 

discussed  previously,  it  is  noted  that  the  basic  nature  does  not 

change  due  to  this  difference  in  initial  conditions.  Also  the  eone 

enthalpy  histories  of  Fig.  30  show  very  little  difftrence  between 

the  lOOK  feet  and  15OK  feet  conditions  and  the  (C_  )  variation  for 

i 

these  is  about  2.5.  Thus  it  appears  that  on  an  absolute  basis  the 
error  must  be  fairly  large  before  any  fundamental  difference  in  wake 
behavior  would  occxu*. 

7.  CONCIUSIONS 

1.  The  analysis  of  Ref.  1  can  be  applied  to  slender  bodies  as  well  as 
blunt  bodies. 

2.  Body  shaping  has  only  a  small  effect  on  initial  core  drag  values, 
these  being  nearly  equal  for  the  sphere  and  cone  of  the  same  base  di¬ 
ameter.  The  Reynolds  number  of  the  body  (size  or  pressure)  is  much 
more  significant  in  determining  the  initial  core  drag. 

3.  The  wane  shock  noy  alter  the  shape  of  the  inviscid  enthalpy  profile 
significantly  for  slender  bodies  (especially  at  low  Mach  number),  but 
does  not  alter  the  vaJ:e  growth  to  any  degree.  This  is  because  the 
general  level  of  the  inviscid  enthalpy  is  so  low. 


COKCLUSIOHS  (Cont'd.) 

4.  The  turbulent  front  growths  behind  the  cones  calculated  are  slower 

X  4 

than  for  blunt  bodies  until  -=r  10  , 

JV 

5.  Axial  enthalpies  drop  very  rapidly  as  the  nose  is  sharpened.  Slender 
cone  wake  enthalpies  decrease  an  order  of  magnitude  faster  than  the 
sphere . 

6.  Bodies  of  the  same  drag  coefficient  but  different  shapes  will  have 
different  wake  histories  until  the  asymptotic  downstream  conditions 
are  reached. 

7*  An  increase  in  core  drag  or  altitude  tends  to  decrease  any  difference 
in  wake  behavior  duo  to  shaping. 

8.  Equilibrium  electron  density  decay  follows  the  fast  drop  of  the 
enthalpy  with  nose  sharpening.  It  was  found  for  the  sphere  -12°  cone 
family  of  bodies  that  the  axial  distance  to  decay  to  a  given  electron 
density  was  directly  proportional  to  the  drag  coefficient.  However, 
this  behavior  Is  not  a  unique  function  of  the  drag  coefficient  for 
all  body  shapes. 

9.  Until  significant  inviscid  wake  drag  is  engulfed  the  slow  cone  wake 
growth  con  cause  the  pure  diffusion  of  a  foreign  species  (e.g.,  ab¬ 
lation  products)  with  zero  recombination  deposited  in  the  wake  to 
take  even  longer  for  the  cone  than  the  sphere. 
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TABI£  I.  Tabulation  of  Body  Paraoeters  -  Re-^otn’  Conditions 


TABLE  II*  Tabulatioii  of  Body  Par«(Mt«rs  -  Ballistic  Bangs  Conditions 


Vake  Sehlnd  Blunt  Body  at  Bypersooic  Speeds 
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2.  Body  Geometry  Schematic 


38 


Invlscid  Flow  Configuration  for  Slender 


Cooqputed  and  Experimental  Shook  Shapes  for  Various 
Sphere-Cones  Including  Wake  Shock 
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Figure  7*  Ckaaputed  Bow  and  Wake  Shock  Shapes  for 

Slender  Cones  (M  =  5»9>  7*8>  and  22) 
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Figure  8.  ConiJarison  of  Theoretical  and  Estperimental 
Shock  Shape  (^/p  =  l/S) 


Figure  9,  Variation  of  n^,  Y.  andfY.  )  vs.  Cone  Angle 
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Figure  10.  Variation  of  n^  and  ve.  Drag  Coefficient 


Figure  13.  lacoa^pressible  and  Compressible  Eoundary  layer  Profiles 
Sefore  and  After  Ease  Region  Expansion 
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Figure  17.  Surface  aud  Wake  Pressure  Distributions  for 
Various  Cones  (M^=  22,  Irt  =  100,000  ft.) 
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Figure  20.  Effect  of  Wate  Shock  on  Enthalpy  Distribution  for 
Conical  Body  (H^=  5*9^  =  •073.) 
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Figure  21.  Effect  of  Wake  Shock  on  Enthalpy  Distrlhufcion  for 
f*  Conical  Body  (M^=  22,  ht  =  100,000  ft.) 
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SLATTERY  8  CLAY  (LINCOLN  LAB)  (REF  4) 


CompaxlBon  of  Theory  and  Experiment  on  Turh'jlent 


Turbulent  Wake  for  a  Blunt  Nosed  and  a  SharpfNosed 
,  Body  With  the  Same 


Figure  89.  Eiuilibri\am  Enthalpy  and  Velocity  Distributions  Along  Walce 
„  Axis  and  Turbulent  Front  (M  =  28,  ht  =  100,000  ft. ) 


Figure  30.  Djuilibrium  Enthalpy  and  Velocity  Distrihutions  Along  Wake 
for  a  Blunt  Hosed  ajai  a  Slirp  Hosed  Body  with  the  same 
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